Quantum Coulomb Blockade in Chaotic Systems 
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We study the Coulomb blockade in a chaotic cavity connected to a lead by a perfectly transmitting 
quantum channel. In contrast to the previous theories, we show that the quantum fluctuations of 
charge, resulting from the perfect transmission, do not destroy the Coulomb blockade completely. 
The oscillatory dependence of all the observable characteristics on the gate voltage persists, its period 
is not affected; however, phases of the oscillations are random, reflecting the chaotic electron motion 
in the cavity. Because of this randomness, the Coulomb blockade shows up not in the averages but in 
the correlation functions of the fluctuating observables {e.g., capacitance or tunneling conductance). 

PACS numbers; 73.23.-b, 73.23.Hk 
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The effect of Coulomb blockade in chaotic systems sets 
up an ideal stage for studying the interplay between the 
quantum chaos and interaction phenomena in a many- 
electron system. By tuning the connection between the 
leads and chaotic cavity. Fig. ^ one can study a rich va- 
riety of non-trivial effects. In the weak tunneling limit, 
discrete charging of the cavity results in a system of sharp 
conductance peaks, which carry information about the 
chaotic motion of non-interacting electrons confined in- 
side an almost closed cavity. In the opposite limit of 
wide channels, charge quantization does not occur, and 
quantum chaos of free electrons in an open billiard may 
be studied. In a broad intermediate region, the charge 
quantization is gradually destroyed, and the chaotic elec- 
tron motion is affected by fluctuations of charge of the 
cavity. The modern experimental technique allows 
one to continuously traverse between these regimes. 

The effect of charging is described by means of the 
Hamiltonian iSl 



virtual transitions of an electron via excited states of the 
cavity Q . The statistics of conductance in this case was 



recently obtained in Ref. |^. 
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where the conventional parameter M is related to the gate 
voltage Vg and gate capacitance Cg hy M — Vg/eCg, and 
Q is the cavity charge. Charging energy Ec for large 
cavities is much larger than the mean level spacing A. If 
the connection of the cavity with the leads is weak, the 
charge is well quantized for almost all M except narrow 
vicinities of the charge degeneracy points (half-integer 
N). The behavior of the differential capacitance of the 
cavity, d{Q) / dVg and of the conductance through the cav- 
ity is quite different for the system tuned to the immedi- 
ate vicinity of charge degeneracy points (Coulomb block- 
ade peaks), or away from these points (Coulomb block- 
ade valleys) . The statistics of the peaks can be related to 
the properties of single electron energies and wave func- 
tions [Q , so that the distribution functions for observable 
quantities can be extracted from the well known random 
matrices results. The transport in the valleys occurs by 
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FIG. 1. Schematic view of a chaotic cavity connected to 
a lead by an ideal quantum channel. Second lead (triangle) 
can be attached to the cavity by weak tunnel junction for 
measurements of two-terminal conductance. 

All the aforementioned results are obtained neglecting 
quatum fluctuations of the charge of the cavity. These 
fluctuations become large with increase of the coupling 
between the cavity and the lead. Then, the difference 
between the peaks and valleys becomes less pronounced 
and eventually instead of the peak structure, one observes 
only a weak periodic modulation. This situation can be 
described neither by the properites of the single-electron 
wave function nor by the lowest order virtual transitions 
via the excited states. The case of almost perfect trans- 
mission of a one-channel point contact connecting a cav- 
ity with the lead was analysed recently by Matveev [Q. 
According to Ref. 0, the Coulomb blockade disappears 
completely if the transmission coefficient of the point con- 
tact is exactly unity and A = 0. The strength of the the- 
ory 0] is in non-perturbative treatment of the Coulomb 
interaction; the drawback is in neglecting the chaotic mo- 
tion of electrons in the cavity. 

In this Letter, we account for both the strong quan- 
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turn charge fluctuations, and the chaotic electron mo- 
tion within the cavity. As we will see, backscatter- 
ing of electrons from the walls of the cavity into the 
channel connecting it to the lead results in oscillations 
of observables with the gate voltage. In the limit of 
perfect channel transmission, the relative magnitude of 
the differential capacitance oscillations is yjK/Ec and 
(A/Ec) In^ {Ec/ A) for the spinless and for spin 1/2 cases 
respectively. If the second lead is attached to the cavity 
(see Fig. |l|) by a weak tunnel junction with conductance 
Go <IC e^/?i, the two-terminal conductance G can be mea- 
sured. The average value of the conductance is supressed 
by the Coulomb blockade, (G) ~ GqA/E'c fluctuations 
of the conductance are of the order of its average. This 
resembles the result for the elastic cotunneling in weak 
coupling regime (5|J^. However, the distinction between 
valleys and sharp peaks is lost. 

Spinless electrons, qualitative discussion - Here, we use 
an analogy to the Fermi liquid theory of the Kondo prob- 
lem due to Nozieres ^ . The dynamics of the system at 
energies smaller than Ec can be well described by the 
lead and the cavity effectively separated from each other. 
Indeed, an electron inside the cavity is not affected by 
the interaction, unless, in the course of chaotic motion, 
it encounters the channel entrance. Even if the trajec- 
tory reaches the channel entrance (Fig. ||), the electron 
is scattered back, provided its energy e (measured from 
the Fermi level) is much smaller than the charging en- 
ergy Ec- The phase shift 5 of this scattering at e = is 
fixed by the Friedel sum rule 6(Q) — irAf, where we used 
the fact that for an ideal transmitting channel the charge 
of the cavity is (Q) = eN 0. On the other hand, the 
variation of the phase of the electron reflection from the 
channel entrance a{e/Ec) = (5(e, A/") — ttM is of the order 
of unity in the energy interval | e | < Ec ■ 

Because an electron at low energy e can not leave the 
cavity, its motion in the cavity is quantized. The one- 
electron energies are periodic functions of J\f, with char- 
acteristic amplitude A and random phases (the latter as- 
sumption will be further elaborated on). The oscillating 
with TV contribution to the ground state energy 5E{N) 
is determined by sum of the eigenvalues in the energy 
interval |e| < Ec\ therefore 5E{J\f) is the sum of approx- 
imately Eel A. ^ 1 random numbers. As the result, the 
average {{6E{N'))) vanishes, and the correlation function 
of the differential capacitance ( in units of Cg) 

Kc{M -M)^ l^dl,dj^{{5E{N)5E{M))) (2) 

takes the form 



In order to explain the functional form of Kc in Eq. (|^) 
and make our argumentation more precise, it is very in- 
structive to evaluate the shift of the energy levels starting 
from the Gutzwiller trace formula |P . The energy of the 
ground state is given hy E = — deN{e), where e is 
measured from the Fermi level, and the integrated den- 
sity of states A''(e) = 9{e — e-i) is expressed as a sum 
over the classical periodic orbits: 



iV(e) = Re^i?j(e)exp 



-Sj{t)^1injb{t,M) 



(4) 



Here Rj is the weight associated with j-th orbit, Sj is 
the reduced action for this orbit. The last term in the 
exponent in Eq. (^) characterizes the reflection from the 
entrance of the cavity, and Uj is the number of such re- 
flections for j'th orbit. We omitted the mean value of 
which is independent of the phase shift (5. In the 
double sum over the periodic orbits, arising in evaluation 
of (^), one can retain only diagonal terms [|lO| because 
different orbits have different actions; the non-diagonal 
terms oscillate strongly and vanish upon the averaging. 
We obtain from Eqs. (||) and (^ 



i^c(A/--)-Re/"j-^El^. 
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where we used the expansion S — S j{e-2) = (ei — £2)^^, 
with tj being the period of j-th orbit. 

In a chaotic system the orbits proliferate exponen- 
tially with the increasing period and cover all the energy 
shell. This exponential proliferation is compensated by 
the damping of the weight which eventually leads 

to the classical sum rule 



dt 

T 



(6) 



valid for periods t much larger than the period of the 
shortest orbit Tl/Et- Energy Et associated with the time 
scale at which the classical dynamics becomes ergodic is 
the counterpart of the Thouless energy for the diffusive 
system. Typically Et ^ Ec , therefore we adopt approxi- 
mation Et Ec A. Apphcation of Eq. (^ to Eq. (||) 
yields 



Kc{f^-) 



'dt^fEct 

t V h 



n^Pr,{t)cos2TTnM_. (7) 



Kc{J\f-) = A(0) ( — ) cos27r7\/l 



(3) 



where the coefficient A(0), see Eq. (^), is the result of the 
effective action theory outlined later. 



Here F (^ 



tt+2ia{e/Ec) 



is a dimension- 
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less function. The number of long orbits is exponentially 
large and they uniformly cover the available phase space, 
therefore, we can employ the stastical description of the 
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reflections by the entrance of the cavity. The proba- 
bihty for an orbit of period t ^ h/ A to reach the en- 
trance even once is small, Pi{T) — tA/h, and one can 
neglect the orbits encountering such a reflection more 
than once. With this simplification, we find Eq. with 
A(0) = dxF{x). The dimensionless integral here con- 
verges at X ~ 1, which corresponds to the orbits with 
t ~ h/Ec ^ fi/ A and justifies the neglection with the 
terms n ^ 1 in Eq. (^. The numerical coefficient given 
by this integral depends on the phase shift a{e/Ec), and 
cannot be found within the simple consideration. 

Finally, we discuss the correlation of the differential ca- 
pacitances as a function of magnetic field. The magnetic 
flux threading a periodic orbit adds a phase (p = AjH/ $o 
to the action in each term j in the Gutzwiller formula (|^) , 
where Aj is the directed area under the trajectory, H is 
the applied magnetic field, and $o is the flux quantum. 
Correspondingly, each term in the sum (^) acquires factor 
cos{Bi Aj /^o) cos{B2 Aj /^o). In a chaotic system, Aj is 
a random quantity which fluctuations grow linearly with 
the period of the orbit, ~ {ETt/h)A^, where 

^ is a shape-dependent area proportional to the area of 
the cavity. Hence, the contributions from the trajecto- 
ries with large period, {{A^j{t))) > ^1/{Bi±B2Y become 
exponentially suppressed. On the other hand, the main 
contribution to the capacitance fluctuations comes from 
the orbits with the period t < fi/Ec- Thus, the correla- 
tion magnetic field Be is controlled by the energy Ec'- 



degrees of freedom of the cavity and obtain a formally 
exact expression for the partition function 



(11) 



where the TV-independent energy fip does not contribute 
to the differential capacitance and can be neglected, and 
Tt stands for the imaginary time ordering. The averag- 
ing (. . .) is performed over the effective one-dimensional 
Hamiltonian 

Ho{N) = ivF I dx {V'I^.Vl - ^Ij^^Vi?} (12) 
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where vp is the Fermi velocity in the channel, and : . . . : 
stands for the normal ordering. The effective action S is 



5= / dTldT2i(ri-r2)^(Tl)V'(T2), (13) 

Jo 

where iP{t) = e'^"° = 0) + iPr{x = 0)) e^''^" are 

the fermionic operators in the Matsubara representation, 
iP{t) = iPH-t). The kernel L(t) can be expressed in 
terms of the exact Green function of the closed cavity: 



L{r) 



Am? 



dydy'^{y)4>{y')dl,,g{T-vy), (14) 
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The exact formula for the correlation function 
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Kc{U;B^,B2) - 
A(z) 
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(8) 

cos 27rA/', 
(9) 



and numerical coefficient in Eq. requires a rigorous 
theory which we outline below. In Eq. (^), we introduced 
the short hand notation B± = Bi ± B2- 

Effective action theory - We start from the Hamilto- 
nian of the system 



H ^ Hp 



He, 



(10) 



where Hp is the Hamiltonian of non-interacting electrons 
in the cavity and in the lead and the Coulomb blockade 
type interaction He is described by Eq. (|I]). 

We calculate the thermodynamic potential il{JV) = 
-TlnZ, Z = Tre"^^, where T = 1/(3 Ec is the tem- 
perature of the system. The second derivative d'j^Vt{N) 
gives the differential capacitance of the system. Because 
the interaction He does not affect the electron dynam- 
ics inside the cavity, we can integrate out the fermionic 



where Q ^ Q — Q^o, the exact Matsubara Green func- 
tion Q is defined with zero boundary conditions (closed 
cavity), Goo is the Green function corresponding to an 
infinite ballistic cavity, and (p{y) is the wave function of 
the transverse motion in the contact. Physically, kernel 
( |T^ ) accounts for all possible returns of an electron to the 
contact during its chaotic motion within the cavity. 

If conductance of the junction connecting the cavity to 
the second lead is small, Go ^ e'^/nh, its effect can be 
accounted for in the second order perturbation theory. 
Derivation similar to that of Eqs. (p^-(p^) yields for the 
two-terminal conductance 



G~o 



-Im 



lim 



*V+iO 



dr- 



sinh ttTt 



n{r) 



, (15) 



where v is the averaged density of states per unit area, 
fin = ZttTtt, is the bosonic Matsubara frequency, and 
function H(t) is given by 



n(^) = 7^7^ f dT,dT2L,{n)L:{T2) 

X {Tre-^F{T)FiO)iiT - Ti)^(r2)). 



(16) 



In Eq. (pq), we retained only the contribution non- 
vanishing in the limit T 0. Unitary operator F in 



Eq. (16) shifts the number of electrons in the cavity by 
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one 11, F''Ho{Af)F = i/o(A/'+ 1); the kernel Le(r) de- 
scribes the motion of an electron from the tunneling junc- 
tion Vt to the entrance of the cavity, and is related to the 
exact one-electron Green function of the closed cavity: 



2m J 



(17) 



In the case of electrons with spins, we imply summation 
over the spin indices in Eqs. (|l|), (H) andjl|). 

The advantage of the representation (|13|) and ( [l^ ) 
is that the Hamiltonian (|2[) is exactly solvable by 
the bosonization technique |M| , and the action S can 
be treated perturbatively provided that A ^ Ec,Et- 
Green function ^ is a random quantity, and its statis- 
tics can be obtained using the well-known expansion in 
the diffuson and Cooperon modes iQ. In the regime 
Ec <C Et the results become universal, i.e. independent 
on the details of chaotic dynamics in the system. Below, 
we present the results of the calculation; details will be 
reported elsewhere p^ ]. 

Spinless electrons, results - For spinless electrons (Zee- 
man splitting exceeds the Fermi energy of electrons in the 
channel) , it suffices to consider only the first order pertur- 
bation theory in S for the calculation of the capacitance 
which gives Eq. (||). Conductance is found from Eq. ( |l5| ) 
and Eq. ( |l6| ) neglecting action S at all. For the average 
conductance, we find 



((G)) = Go 



A 



^Ec 



A(0), 



(18) 



where 7 « 1.78... is the Euler constant and function 
A(a;)is defined in Eq. The correlation function of the 
conductance fluctuations is given by 



{{SG1SG2)) 
((G))2 



Bl 



cos^ ttTV- 



A2(0) 



(19) 



where the correlation magnetic field is given by Eq. (|^), 
we introduced the short hand notation, G,; = G{J\fi,Bi), 
and B± = Bi ± B2, A/l = A/l - A^2- 

Spin 1/2 case, results - In the spin 1/2 case, the non- 
vanishing correlation function of the differential capaci- 
tance appears only in the second order of perturbation 
theory in S. In zero magnetic field, we find 
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Stt^E^ 



In* 



Ec 
T 



cos 27r7V- . 



(20) 



We present here also the correlation function in the uni- 
tary limit (5+ ^ Be), 
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Br 
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Ec 
T 



B^ 

s7 



cos27r7V-, (21) 



valid for the region of fields Bc^T/Ec < i?- < i?c- 

To calculate the average conductance, one can neglect 
S, similarly to the spinless case. This yields 



(22) 



The main conribution to the conductance correlation 
function is independent on M- , but depends on the mag- 
netic field; in the unitary limit 



In 



(23) 



The characteristic magnitude of oscillating with TV con- 
tribution to the correlation function is smaller than 
Eq. ( p3| ) by a parameter A/£'c. 

At low temperatures T, the above results diverge, 
which indicates that the higher order terms in S should 
be taken into account. Using a variational approach [p^ , 
we conclude that T should be replaced with A In^ {Ec /A) 
in Eqs. (|ol)-(|3l). 

In conclusion, Coulomb blockade oscillations persist 
even if a cavity is connected to a lead by a perfect single- 
mode channel. This effect is ignored if one approximates 
the electron spectrum in the cavity by a continuum [Q, 
or neglects the charge quantization [|l5| . 
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